Mon. Not. R. Astron. Soc. 0Q0.[T1[m(2Q05) Printed 5 February 2008 (MN LffeX style file v2.2) 



Voigt Profile Fitting to Quasar Absorption Lines: An Analytic 
Approximation to the Voigt-Hjerting Function 



Thorsten Tepper Garcia* 



\0 

o 
o 

<: 

0^ 



> 

(N 

(N 
O 

o 

O 



Institut fiir Astrophysik, Georg-August Universitdt, Friedrich-Himd-Platz 1, D-37077 Gottingen 



Accepted — . Received — ; in original form — . 



ABSTRACT 

The Voigt-Hjerting function is fundamental in order to correctly model the profiles of ab- 
sorption lines imprinted in the spectra of bright background sources by intervening absorbing 
systems. In this work we present a simple analytic approximation to this function in the con- 
text of absorption line profiles of intergalactic HI absorbers. Using basic calculus tools, we 
derive an analytic expression for the Voigt-Hjerting function that contains only fourth order 
polynomial and Gaussian functions. In connection with the absorption coefficient of inter- 
galactic neutral hydrogen, this approximation is suitable for modeling Voigt profiles with an 
accuracy of 10 '* or better for an arbitrary wavelength baseline, for column densities up to 
Nhi = 10^^ cm"^, and for damping parameters a < 10"^, i.e. the entire range of parameters 
characteristic to all Lyman transitions arising in a variety of HI absorbing systems such as 
Lya Forest clouds, Lyman Limit systems and Damped Lyo- systems. We hence present an 
approximation to the Voigt-Hjerting function that is both accurate and flexible to implement 
in various types of programming languages and machines, and with which Voigt profiles can 
be calculated in a reliable and very simple manner 

Key words: methods: analytic, quasars: absorption lines, line: formation, line: profiles, line: 
identification 



1 INTRODUCTION 

Absorption processes and their signatures (absorption lines) im- 
printed on the spectra of bright background sources (quasars. 
Gamma-ray bursts, etc.) are one of the main sources of information 
about the physical and chemical properties of intervening systems. 
It is well known that information about their temperature, density, 
chemical abundances, and kinematics can be extracted from the 
analysis of these absorption lines. For instance, a detailed insight 
into the physical state of the intergalactic medium (IGM) is pro- 
vided by the analysis of the absorptio n lines found in the spec- 
tra of distant quas ars (QSOs) (see e.g. IHu et a/.lll9 95: Kim et al. 
11997 *. 2001, 2002"). These lines are due mainly to absorption by 
neutral hydrogen (HI) present in a class of low column density 
absorbers generally known as Lyor Forest, and due to other ele- 
ments in low ionisation stages (CII, CIV, Sill, Mgll, Fell, Oil, etc.), 
which arise in higher column densities absorbing systems associ- 
ated with galaxies, such as the Lyman Limit Systems (LLSs) and 
Damped Lyo- Absorbers (DLAs). A wealth of information about 
the distribution, density, temperature, metal content, etc. of these 
systems is now available as a result of exhaustive and extensive 
studies of QSO absorption lines (see e.g. ■Rauch.l998:.Raa.2005t 
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for a review on Lyo- absorbers, and on metal systems and DLAs, 
respectively). 

In this type of analysis, and within the realm of a given cos- 
mological model, the number and observed central wavelength of 
the absorption lines provide information on the spatial distribution 
of the absorbing systems. Furthermore, knowledge about the actual 
physical state of these systems can be obtained basically from the 
line profiles. Both line counting and line profile measurement are 
tricky tasks though, since the accuracy with which they can be per- 
formed highly depends on the resolution of the observed spectra. 
For instance, depending on the spatial distribution of the absorbing 
systems, lines can appear very close to each other or even super- 
pose (line blending), and a low spectral resolution may lead to the 
misidentification of the resulting composite profile as being a sin- 
gle, complex one. Because of this same reason, the determination 
of the exact shape of each individual absorption profile is far from 
being trivial, and misidentified profiles may lead to wrong conclu- 
sions about the properties of the absorbing systems. 

If one assumes that the physical state of the absorbing medium 
is uniquely defined by its temperature and column density, sin- 
gle absorption line profiles are ideally described by Voigt profiles. 
Mathematically, a Voigt profile is given in terms of the convolu- 
tion of a Gaussian and a Lorentzian distribution function, known 
as Voigt-Hjerting function iHiertine 1938), and a constant factor 
that contains information about the relevant physical properties of 
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the absorbing medium (cf. Sect. 12. 1> . Any departure from a pure 
Voigt profile in the observed hues is expected to yield information 
about the kinematic properties (non-thermal broadening, rotational 
or turbulent macroscopic motions), as well as spatial information 
(clustering) of the absorbing systems. 

The Voigt-Hjerting function has long been known and, con- 
sequently, various numerical methods to estimat e and tabulate 
this function have been developed and presented (Hierting 1938; 
[Har ris 1948; Finn & Muggles tone 1965). Also, a great effort has 
been done in order to derive a semi-analytic approximation to this 
function that reproduces its behavior with high accuracy (see e.g. 
[whiting 1968; KielkoDf 1973; Monaghan 1971), being the latter by 
far the one with the highest accuracy. With the help of these meth- 
ods, computational subroutines have been developed that make it 
nowadays possible to nu merically integrat e this function for a wide 
parameter space (see e. g. iHumlfcekl 1983) . 

The aim of this work is to make a further contribution to the 
practical handling of the Voigt-Hjerting function in order to com- 
pute Voigt profiles. Starting with an exact expression for this func- 
tion in terms of Harris' infinite series, we argue why this series may 
be truncated to first order in a in the context of intergalactic HI ab- 
sorption lines. We then show that the second term of this series can 
be approximated with a non-algebraic polynomial function, which 
is mathematically simple to handle in the sense that is does not 
contain singularities. Such an analytic, 'well-behaved' expression 
in terms of simple functions as presented here is very attractive, 
since it allows one to replace the many steps and operations needed 
for numerical integration, or to read from look-up tables of values, 
by a single line with simple operations. It is also extremely flexible 
to implement in various types of codes and machines, and is par- 
ticularly useful for computational routines in higher-level program- 
ming languages {e.g. IDL, Mathematica, Maple, etc.), in which nu- 
merical integration or look-up table reading is cumbersome, espe- 
cially if absorption line profiles have to be calculated many times 
with moderate precision and relative high speed. For instance, such 
an analytic expression should be very useful to synthesise Lya ab - 
sorption spectra as in e.g. .Zhang et a/.,.l.l997.) ; Richter et a/.li2005h . 
or in line-fitt ing algorithms like AUTO VP jpave ef a/.lll997l) or 
FITLYMAN ('Fonta na & Ballestej[T993) . used to obtain line pa- 
rameters such as redshift, column density, and Doppler width from 
absorption Voigt profiles imprinted on observed spectra. 

In the next section we briefly outline the origin of the Voigt- 
Hjerting function in the Physics of absorption processes, and define 
the context in which the desired approximation of this function is 
of interest to us. In Sections|^and|31we derive this approximation, 
and in Section|5|we compare the accuracy and speed of a numerical 
method for computing Voigt profiles based on our approximation to 
other existing methods. In Section |S| we present an application of 
our method to model Voigt profiles, and in Section0we summarise 
our main results. 



2 THE VOIGT-HJERTING FUNCTION IN THE 
CONTEXT OF HI ABSORPTION LINES 

2.1 The Absorption Coefflcient 

The probability of a photon with an energy E = he/ A to he ab- 
sorbed within a gas with column density N and kinetic temperature 
T, also known as absorption coefficient, is given by 

Ti(A) = (Ci-'N-a)-H(a,x(A)), (1) 

where 



4V^^ 
m^c 1 ,- 

is a constant for the ;th electronic transition caused by the pho- 
ton absorption. Here in^ is the electron mass, / is the oscillator 
strength, and F, the damping constant or reciprocal of the mean 
lifetime of the transition. The function H is the so-called Voigt- 
Hjerting function and is given by 

a r+" e->' 

H{a,x)=- 7j— ^dy. (2) 

Let Aj = h c/Ej be the resonant wavelength of the corresponding 
transition , and AAo = ^/l; the thermal or Doppler broadening, 
which defines a Doppler unit. Here, the Doppler parameter b is re- 
lated to the kinetic temperature of the gas via b = -^2kT/mp, where 
k is the Boltzmann constant and nip is the proton mass. It follows 
from these definitions that the damping parameter 

A}r, 

a = 

4;rcA/lo 

quantifies the relative strength of damping broadening to thermal 
broadening, and that the variables x = and y = I are just 

the wavelength difference relative to the resonant wavelength in 
Doppler units, and the particle velocity in units of the Doppler pa- 
rameter, respectively. 

The particular form of the absorption coefficient Q induces a 
characteristic absorption feature known as Voigt profile. Hence, the 
Voigt profile, and consequently the Voigt-Hjerting function, nat- 
urally arise in the process of absorption-line formation, when one 
assumes that the physical state of the absorbing medium is uniquely 
defined by its density and kinetic temperature. Generally speaking, 
it is the physical conditions what determines the shape of the ab- 
sorbing features, i.e. the line profiles. Conversely, it is true that line 
profiles give information about the physical state of the absorbing 
medium. In particular, line profiles other than Voigt profiles give in- 
formation about the departure of the physical conditions assumed 
here. 

The class of HI absorbers present in the IGM (Lycf Forest 
Clouds) and associated with galaxies and larger structures (LLSs, 
DLAs) can be characterised by their column density Nhi and ki- 
netic temperature, and consequently their absorption features ob- 
served on e.g. QSO spectra are well described by Voigt profiles. 
Their observed column densities span a range of ten orders of mag- 
nitude, approximately from 10'--10-^ cm"^, and have temperatures 
that correspond to Doppler parameters approximately in the range 
10 - 100 km s"', with a median value around b„, = 36kms"' that 
decreases with redshift I Kim et al. 1997). For such a range in b, the 
damping parameter a for the Lyman transitions of intergalactic HI 
spans a range of 9.3 ■ 10"' - 6.05 • lO""*. In this case, high values 
of a correspond to Lya, while lower values are typical for higher 
order Lyman transitions. Transitions of other elements, such as C, 
Si, Mg, Fe, O, etc., and their various ionisation stages, are also typ- 
ically found on QSO spectra, and their damping parameters cover a 
range which strongly overlaps with that of HI. This can be seen in 
FigureQ where we show the distribution of a for different elements 
(including HI) in different ionisation stages, for a Doppler parame- 
ter ^Hi = 36 km s"' for HI, and assuming that the doppler parameter 
for other elements is related to bni via bx = bm ■ V™hi/"^x> where 
ntx is the mass of element X. For clarity, we have grouped all dif- 
ferent ionisation stages of a given element under a single label. The 
values of the atomic constants (central wavelength of the transition 
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Figure 1. Value of the damping parameter a, assuming bm = 36kms"', 
for different elements in different ionisation stages typically found in QSO 
spectra. The doppler parameter for other elements is assumed to be given 
by bx = bm ■ y/ mm I nix , where mx is the mass of element X. Elements are 
listed by increasing atomic mass on the x-axis, and the logarithmic value 
of a is given on the y-axis. For clarity, different ionisation stages of a given 
element have been grouped under a single label. The shaded area marks the 
range above the largest value of a for the intergalactic HI Lyman transitions. 
Note that for the sake of completeness, other elements than those found to 
date in QSO spectra have been included as well. 



and F- value) have been taken from iMortonI i2003h '. The shaded 
area contains the values of a above the range characteristic to in- 
tergalactic HI for which our approximation to the Voigt-Hjerting 
function-derived in the next sections-cannot be applied or should 
be applied with caution. Note, however, that the region spanned by 
a for intergalactic HI, i.e. the region underneath the shaded area, 
contains most of the a-ranges spanned by all other elements, espe- 
cially those corresponding to C, O, Mg, and Fe. For other elements, 
such as Cr or Zn, the damping parameter has values right at the up- 
per limit of this range. For this reason and for the sake of simplicity, 
in the following we will constrain our discussion to the Lyman ab- 
sorption lines of intergalactic HI, but the reader shall bear in mind 
that the the discussion and method to synthesise Voigt profiles pre- 
sented in this work can be directly applied to the transitions of other 
elements associated with HI absorbing systems, according to Fig- 
ure Q 



2.2 The Absorption Coefflcient of HI at Low Column 
Densities 



Following iHarrigj il948!), it is true that for a < 1, i.e. when 
Doppler broadening dominates over damping broadening, the 
Voigt-Hjerting function can be expressed as 



H{a, x) = Yj (^) a" , 

where the functions /f„(x) are defined by 

H„ (x) = -^^^ r v"e-<"'2'' cos(xv) dv . 
wn! Jo 



(3) 



(4) 



' A list containing the values of the damping parameters for 
the elements and their different ionisation stages shown in Fig- 
ure Q] is available in plain-text format at www.astro.physik.uni- 
goettingen.de/~tepper/hjerting/damping.dat. Please consult this list in order 
to know the exact value of a for a given element in a given ionisation stage. 




Figure 2. Departure of the Voigt-Hjerting function from a pure Gaussian 
function for a = 10"^ (solid line) and a = IfT^ (dashed line) as a func- 
tion of X. The departure is given as the logarithmic difference between 
the Voigt-Hjerting function and the zeroth order approximation in the form 
SHq = \ - HqI H . The greater this quantity, the less accurate is the zeroth 
order approximation. 



These functions are bounded with respect to n and x, and they have 
values of the order of unity. Indeed, taking the absolute value of the 
integral, neglecting the cosine, which takes values of the order of 
unity, and computing the resulting integral one can show that 



\H„{x] 



^ — « 1.123, 

-0r 



(5) 



for all n 6 {0, 1, 2, . . .) and x 6 R. From this it follows that if a <c 
1, the Voigt-Hjerting function can very well be approximated to 
zeroth order in a by the first term of the series ^3}' H(a, x) x 
Ho(x), as first noted bv IWalshawl h955h . Note that this result is 
exact in the limit a ^ 0. Taking the definition 0, it follows that 



Ho(x) 



(6) 



and thus H(a, x) « e"'", for a <k 1. We call this the Voigt-Hjerting 
function to zeroth order. 

But what actually means that the condition a <k 1 be satisfied, 
so that this zeroth order approximation can be safely used to model 
absorption line profiles? In order to address this, we compute for 
the extreme values of a for intergalactic HI, a a 10"^ and a a 10"*, 
the departure of the Voigt-Hjerting function from a pure Gaussian 
function. This is shown in Figure |5| as the logarithmic difference 
between^ H and Hq relative to // as a function of x, i.e. the quantity 
6Ho = I - Ho/H. Note how the zeroth order approximation com- 
pletelly differs from the actual Voigt-Hjerting function at x > 3.5 
for a = 10"'', and at x > 4.5 for a = 10"*. Thus, even a value of 
a as small as 10"** is not a sufficient condition for the zeroth order 
term to be a good approximation of H for an arbitrary range in x. 

In the context of the HI absorption coefficient, the condition 
on a for which the approximation to zeroth order of the Voigt- 
Hjerting function is valid, translates into a restriction on Nhi, and 
more specifically on the quantity C, • a • Nhi. This can be seen by 
taking a glance at equation 0: Despite the fact that the condition 
a <s: I holds, the product C/ ■ a ■ Nhi can be very large for high 



^ The values of the Voigt-Hjerting function were computed numerically us- 
ing a routine based on Monaghan's algorithm in ^Murphv 1 2QQZ) (cf. Section 
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enough'' column densities, and since the functions H„ are of the or- 
der of unity, such terms can have a significant contribution to the 
absorption coefficient. Along this line of reasoning, and taking into 
consideration that the constant C, in equation Q is of the order of 
10~" cm^ for all Lyman transitions, it is clear that the departure of 
the Voigt-Hjerting function from its zeroth order approximation be- 
comes significant in the wavelength ranges x > 3.5 and x > 4.5 (in 
Doppler units) at column densities Nhi > 10'^ cm"^ for a damping 
parameter a = lO"*, and at column densities Nhi ^ 10" cm"^ for 
a = 10"**. Since intervening HI absorbers typically have column 
densities in the range 10'^ - 10^" cm"', a Gaussian approximation 
to H for modeling intergalactic HI absorption line profiles is only 
suitable for the low end of the column density range. 

In addition to the factor C, • a • Nhi being large and even 
more decisive for the uselessness of the zeroth order approxima- 
tion for an arbitrary range in x, is the fact that the zeroth order term, 

2 

, rapidly decreases for large values of x and is therefore over- 
whelmed by higher order terms in Harris' expansion, which hence 
dominate the behavior of the absorption coefficient, as already seen 
in Figure|2| To shed some light on this, consider the following nu- 
merical example: Out to x s; 4 (in Doppler units), and for a = 10"'', 
the Voigt-Hjerting function is of the order of 4.2 • 10"*. The ze- 
roth order term in the series J3} satisfies Ho (x = 4) x 1.1 • 10"', 
whereas the first order term (a ■ Hi) (x = 4) » 3.9 ■ 10"*. Thus, 
at large enough x, the behaviour of t is evidently governed by the 
terms of order I'm the series jsj- 

2.3 Higher Column Densities and First Order Term 

Due to the arguments stated above, and even though the damping 
parameter satisfies a <c 1 , it is clear that the absorption coefficient 
of intergalactic HI cannot simply be approximated by a constant 
times e"'' for an arbitrary wavelength range and column densities 
grater than lO" cm"^. One actually has to take into account terms of 
higher order in the series at least to first order in a, i.e. H(a, x) a: 
(Ho + a- Hi) (x) . In fact, one should take into account all terms up to 
Mh order for values of (C, • a ■ Nhi)" ' that are nearly equal or greater 
than the absolute difference between the sum 2lI=o (-'■) ' and 
the exact Voigt-Hjerting function. However, as we shall show next, 
the approximation to first order in a in Harris' expansion is enough 
to model absorption line profiles with moderate to high accuracy for 
the range of parameters (a, Nhi, C,) characteristic to intervening HI 
absorbers. 

In order to prove the above statement, we look at the relative 
contribution of the zeroth and first order terms to the Voigt-Hjerting 
function for a e {10"**, lO"**) and an arbitrary range in x. This is 
achieved, for example, by computing the logarithm of the quantity 
6Hi = |1 - (Hq + a ■ Hi)/H\ as a. function of x for the two extreme 
values of a, as shown in Figure|3| Here we have advanced the func- 
tion Hi which is being handled in the next section. Note that the 
greater the contribution from the zeroth and first order term to H, 
the smaller the quantity SHi. In this case, as can clearly be seen, 
6Hi takes on values of the order of 10"' or less over the whole 
wavelength range shown here and for the whole range in a for in- 
tergalactic HI. 

If one takes into account that H(x) < 1 for all x, it is obvi- 
ous that the relative difference is equal or greater than the abso- 
lute difference, i.e. SHi > \H - {Ho + a ■ H\)\. Thus, in the case of 

' 'High enough' means in this case that the condition Nhi > (C, ■ a)"' is 
satisfied. 
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Figure 3. Contribution from the zeroth and first order terms of the series J5) 
to the absorption coefficient, when (C, ■ a ■ Nhi) is of the order of unity. The 
curves show the logarithmic difference between H and Ho + a ■ H\, relative 
to // as a function of .v for a = 10"^ (solid line) and a = 10"'* (dashed fine). 
The smaffer this dift'erence, the greater the contribution from the zeroth and 
first order term to H. Note that the behavior of H is indeed dominated with 
a difference of seven or more orders of magnitude by the first two terms 
of HaiTis' expansion for the whofe range in x shown here. The vafues for 
Hi were cafcufated according to equation JS), Section |3] using numericaf 
integration to compute the function F. 

a X 10"'* and 5Hi = 10"', the departure of H from its first or- 
der approximation becomes significant at column densities Nhi > 
(a ■ SHi ■ C,)"' = 10^" cm"', and at even larger Nhi for a < 10"* 
and/or smaller SHi. Hence, the first two terms of Harris' expan- 
sion dominate the behaviour of H over the whole wavelength range 
shown here, with an accuracy of 10"' or greater, for the range in a 
characteristic to intergalactic HI. On this basis, we consider thar an 
approximation to first order in a of the Voigt-Hjerting function in 
terms of the functions Ho and Hi is suitable to model Voigt profiles. 
Since the function Ho is known and simple per se, we now turn to 
the task of finding an approximation to the function Hi in terms of 
a simple, analytic expression. 



3 THE DAWSON FUNCTION REVISITED 

According to definition ^4} we have 

-4 f" _„2 

Hi (x) = — — I V e cos(2xv) dv . (7) 
-0r Jo 

Integrating this equation partially, and computing the resulting Si- 
nus transform of a Gaussian it follows iMihalaJl97flh 

//,(x)= ^[l-2xF(x)], (8) 

where we adopt the notation first introduced by iMiller & Gordoi] 
il93lh 

F(x) = e"'' r e'''dv. (9) 
Jo 

F is known as the Dawson function iDawsonll898l) . 

It is evident that finding an approximation for Hi translates 
into the same problem for F. We thus want to show that a simple 
analytic expression can be found which approximates the Dawson 
function, and hence the function Hi , accurately enough, in order for 
an approximation of the Voigt-Hjerting function in terms of these 
functions to be useful for synthesising Voigt profiles. 
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3.1 Properties of the Dawson Function 

We briefly want to state some important properties of the Dawson- 
Function. First, this function is antisymmetric, i.e. F{-x) = -F(x) 
for all X € R. Because of this, from now on we restrict our anal- 
ysis to X > 0. Besides, it has no roots in the positive semi-axis, 
and F(0) = 0, as can easily be seen from its definition (Fundamen- 
tal Theorem of Calculus). Furthermore, F is bounded, since Hi is 
bounded itself (cf. Sect. l2.2t . Indeed, differentiation with respect to 
X gives 



— F{x) = l-2xF{x). 
dx 



(10) 



Hence, the upper bound is given by F(xq) = (2xo) ' , where xq is de- 
fined by equating F's derivative to zero and solving for x. Actually, 
F has it s maximum at xp = 0.92413 with F(0.92413) = 0.54104 
(see e.g. Abramowitz & Steaun eds. 1965). From this last equation 
it follows also that equation (S} can be rewritten as 

Hi(x)= ^-^F{x). (11) 
^Jn dx 

We want to know how the function F behaves asymptotically, 
i.e. for X <c 1 as well as for x » 1. Using the power series of the 
exponential we get 



F(x). 



dv 



n\ In + 1 



(12) 



In this case, the sum and the integral operator commute, since the 
power series of the exponential converge s uniformly in any interval 
[a, b], particularly for v 6 [0, x] (see e.g. Forster 1983). Expressing 
the term e""' by its corresponding power series and rearranging 
terms, this last equation reads in explicit form 

F(x) = x-(l-^x^ + ^x* + ---). (13) 

Thus, for X <K 1 the Dawson function behaves asymptotically up to 
third order as 

2 , 

(14) 



X « 1 . 



F(x)«x-(l--x2), 

In order to investigate how F behaves for x » 1, we first rewrite 
the function J9} with the replacement v' = x - v and the aid of the 
power series of the exponential as 



with the definition 

x)= { v-"e-2" 
Jo 



/„( 



(15) 



(16) 



It is not hard to see that every term of the series <15> is separately 
bounded with respect to x as well as n. Making the replacement 
v' = 2vx in eq. <16> . and integrating for x ^ we get for n 6 No 

1 



/„(x) 



-{2n)\{l e-'' R„(x), 



(17) 



(2x)2"+i 

where R„{x) is a rather cumbersome polynomial function. Now, for 
X s> 1, we may drop all terms which contain an exponential factor 
and in this way we get the asymptotic form 



F{x)^y ^(2n)!, x»l. 

^ ' Zj (2x)2"+i ^ ' 



(18) 



It can be seen from this expression that F vanishes as (2x) ' for 
X — > CO, and that the first derivative <10> . and thus the function Hi, 



also vanish in this limit as (2x^)"'. From equation jlO> it is also 
true that F' converges to unity for x — > and that Hi converges to 
(-2 / ^fH) in this limit. Since we want our approximation to F, and 
consequently to Hi and H, to be valid in the whole range x 6 [0, oo), 
we require it to fulfill both these conditions as well. 



4 THE ANALYTIC APPROXIMATION 

Equation <15> . together with eq. <16t . represent indeed an exact ex- 
pression for the Dawson function. However, these expression are 
not suitable for practical computation. We therefore explore the 
possibility of finding an analytic expression which is easy to handle 
and which can be used to compute the value of F{x) for x 6 [0, oo). 
In particular, we shall see if it is possible to truncate the series <15> 
in order to find an approximation to F, which has all its properties 
(antisymmetry, boundedness, etc.), which converges for x — > oo as 
well as for x — > 0, and which is well defined in the whole range 
[0, oo). For instance, equations <14> and <18> do not fulfill these re- 
quirements. Nevertheless, they show us how our desired function 
has to behave asymptotically. 
Let us define 

Z)„(x)= V -/„(x). (19) 
where the /„'s are given by eq. <16> . Using this definition we get 

Diix) = i,-.-^)\^^.^\-.-^^' (20) 

It is easy to show that this function behaves qualitatively in the 
same way as F does, i.e. it is antisymmetric, bounded, and has no 
roots in the positive semi-axis. Furthermore, both these functions 
have the same asymptotical behavior. Indeed, up to third order we 
have for X <K 1 

2 



Di{x)^ x-{\-^x\ 



X « 1 , 



(21) 



as can be shown by expanding the exponentials in eq. j20> in terms 
of their power series. A glance at eq. <14t makes the similarity be- 
tween F and Di evident in this limit. For x » 1 we get from eq. 
<20> . neglecting the exponentials. 



X » 1 . 



that is the same as for F (see eg. llSt . Furthermore, the first deriva- 
tive of Di with respect to x is unity at x = and vanishes as (2x-)"' 
for X — > oo. However, the maximum of Di is at x = 0.87269 with 
£»i (0.87269) = 0.52212, i.e. at slightly diff'erent values from those 
of F. 

The magnitude and range of the error in approximating the 
Dawson function by the function Di can be estimated, at least qual- 
itatively, in the following way: From equations {15\ and J19> it fol- 
lows that 



F(x) = £),(x)+ V^/„(x). 



(22) 



The sum in this last equation, i.e. the absolute error in our ap- 
proximation F(x) X Diix), is a positive semi-definite quantity, 
since each term of the sum has this property. Hence, it is true that 
< F'i(x) < F(x) for X e [0, oo). Furthermore, since both F 
and Di converge to zero for x — > as well as for x — » oo, and 
both these functions are bounded, it follows from eq. H2\ . that 
the error also vanishes for x ^ as well as for x — > oo, that 
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it is also bounded, and that its maximum value is less equal than 
max{F(x) — Di(x)}^q. From this, and since their respective max- 
ima are slightly shifted with respect to each other, one is led to the 
conclusion, that the error is constrained to a narrow range in x and 
that the maximum error in our approximation to F occurs near the 
maxima of these functions, i.e. in the vicinity of x= 1. 

We wont further try to quantify the actual error in our approxi- 
mation to F. It shall be enough to know, for our purposes of finding 
an approximation to the Voigt-Hjerting function, that the error in 
the approximation F(x) x Di(x) is bounded and constrained to a 
narrow wavelength interval around x = 1. Besides, we will indi- 
rectly estimate the error in this approximation when quantifying 
the error in our approximation to H in Section|5| 



4.1 The Voigt-Hjerting function to First Order 

Once we have found an approximation to the Dawson function, 
we can use it to give the desired expression for the Voigt-Hjerting 
function using Harris' expansion to first order in a. Replacing in 
equation <1 U the function F by our approximation D[ (eq. I20> . 
taking the corresponding derivative, and rearranging terms we get 



Hi(x)x -^K(x)e-'^ , 
where we have defined 



K(x): 



1 

2? 



{Ax^ + 3) (x^ + 1) e^ 



— (2x^ + 3)sinhx2 
x^ 



(23) 



(24) 



We want to highlight the fact that equation <23> is well defined, 
i.e. it has no singularities in the whole interval [0, oo). Further- 
more, it converges to the correct value in the limits x — > and 
X — » oo. Indeed, it is easy to show that limj-_,o //i (x) = ^ , and 
lim,-^„ Hi{x) = 0, whether one uses for H[ the exact expression 
Jsj or the approximation <23> . In contrast, in approximations to the 
Voigt-Hjerting fun ction to model Voigt profiles ofte n used in the 
literature (see e.g. ISoitzeJIigyst IZhang et a/.lll99'il and given in 
the form Ci ■ e"* + C2 ■ \ , where the c,'s are constants, the second 
term which represents the Lorentzian damping clearly diverges for 
X — > 0, and one has to artificially define the wavelength range in 
which this second term is used. It is customary to neglect this term 
for low column densities and in the vicinity of x = 0, but how to 
exactly choose the radius of the vicinity is not clear and completely 
arbitrary. However, with an expression like equation )23t at hand, 
no such assumptions have to be made. 

Taking into account that a <c 1 in order to neglect terms of 
order n > 2 in the series J^, we get, using the expressions for 
Ho (eq.|5| and Hi (eq. l23L that the Voigt-Hjerting function to first 
order in a is given by 



H{a, x) a e 



■ a— K(x) 



(25) 



This expression is symmetric in x, as it should be, and thus it is 
valid for X 6 R and a « 1 . According to this equation, the Voigt- 
Hjerting function can be regarded as a "corrected" Gaussian func- 
tion, where the correction term depends on the parameter a. In the 
context of the absorption coefficient of HI, this correction term also 
depends on the column density Nhi, of course, via the quantity 
fl-Nm. 



5 ANALYSIS 

In order to quantify the quality of our approximation to H, we per- 
form a test on speed as well as on precision, comparing a numer- 
ical method to compute H, based on our approximation, to other 
standard, available methods to numerically compute this function. 
For this purpose, we use the approach and c orresponding computa- 
tional routine developed byl Mur p^ i2002h . which consists of the 
numerical implementation in FORTRAN of four different methods 
to compute H: Harris' Hi and H2, Humlicek's, and Monaghan's. 
In Murphy's notation, Harris' H\ and H2 correspond to the Voigt- 
Hjerting function approximated by the first three and five terms of 
the series expansion J^, respectively. Humlicek's optimized algo- 
rithm and Monaghan's differenti al app roac h to approxi mate H are 
explained in detail in Humlfcek 1 1982) and Mo nagha re- 
spectively. Our method to compute H consists simply in the numer- 
ical implementation in FORTRAN* of equations I24> and <25> . 



5.1 Speed 

Following Murphy's approach, the relative speed of all five 
methods are determined by calculating the time that a rou- 
tine based on each method requires to compute the Voigt- 
Hjerting function for x 6 [0, 10] and damping parameters a in 
the range 10"** - 10"* for a total of 1.5 • 10' runs. In this 
way, we get that the relative speed"" of each method in the or- 
der Hi : (this work) : Humlfcek : HI : Monaghan corresponds to 
1 : 4.2 : 5.8 : 6.8 : 66.6, independent of a. According to this result, 
our method is second fastest. 



5.2 Precision 

The precision of our method relative to the other methods men- 
tioned above is determined in the following way: First, the values 
of the Voigt-Hjerting function computed using Monaghan's algo- 
rithm for X 6 [0, 10] and a e [10"**, 10"*] are taken as fiducial. 
Then, the value of H for the same range in x and a is computed us- 
ing each of these methods, and the logarithmic difference between 
each method and its fiducial value, i.e. the quantity log|Q 6H with 
dH = 1 - H^fiijoj/Huon, is calculated as a function of x for each 
different a. The result is shown in Figure |4] It can be seen from 
this figure that Harris' H2 is the second most precise method to 
compute H, if one takes Monaghan's algorithm as fiducial, but is 
six times slower than Harris' Hi, and one-and-a-half times slower 
than our method, as stated in the previous section. Our method has 
a precision of 10""* or better for x < 4. For all values of a, the dif- 
ference peaks around x = 4 to a value of the order of 0.01, and the 
precision increases again for values of x > 4. The precision is better 
for smaller a, as expected, since the zeroth order term gains in im- 
portance in our approximation for decreasing a. For a < 10"* and 
X < 1.5, our method is more precise than Harris' Hi or Humlicek's, 
and, as seen above, 1.5 times faster than the latter. 



* The rearrangement of equations 1241 and <251 that leads to the smallest 
number of operations reads, in code syntax, 

H(a,.x) = Ho-a I -Jn I]? - [HQ-Ho-iA-j? ■ y? +1 ■ + A + Q) - Q - \\.. 

where the terms x- , Hq = e" and Q = 1 .5 ■ have to be computed just 
once. 

^ This calculations were performed on an Intel Xeon 3.2 GHz processor. 
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Figure 4. Precision of different metliods to compute tlie Voigt-Hjerting function, relative to Monaghan's differential method. Shown here is the logarithmic 
difference as a function of x, i.e. the quantity logig 5H, with 6H = 1 - H^ethodlHMon, where Hmou and Hi„g,i,oj are, respectively, the Voigt-Hjerting function 
computed using Monaghan's algorithm and each of the methods mentioned in the text: HI (solid line), H2 (short-dashed line), Humh'cek (dot-dashed line), 
this work (long-dashed line). Each panel corresponds to a different damping parameter. Here we chose the range of a characteristic to intergalactic HI, i.e. 
a 6 [10^^, 10"'']. This graph was created using the approach and con'esponding routine developed by MuiTihv ( 2002 ), and is adapted from Figure A.l of the 
same work. 



5.3 Modeling of HI absorption profiles 

We now turn to analyse how accurate is our method in order 
to model HI absorption profiles. Taking the whole range in col- 
umn density logNni 6 [12.0,22.0] dex and Doppler parameters 
b e [10.0, 100.0] kms"' characteristic to intergalactic HI, we syn- 
thesise for each pair (Nhi, b) (with a resolution of A log Nhi = 0.05 
dex, and Ab = 0.5 km s"') a single Lyo- absorption profile in the 
range A e [100, 2300] A with a resolution of A^ = 0.01 A. The ab- 
sorption profile is synthesised according to equations Q and J2j, 
using both our method and Monaghan's to compute H. We then 
compute for each pair (Nhi, h) the absolute value of the difference 
between the profiles generated using these two methods relative to 
Monaghan's, i.e. 6V = \l - e^'^^l, with At = t„,„. - t^o,,, as a func- 
tion of wavelength for the whole wavelength range, and pick the 
maximum value of this difference in the entire range. We choose to 



do so in order to pick up the worst cases possible, i.e. those with the 
lowest accuracy, and put in this way a stringent lower limit to the 
accuracy of our method. The result is shown as a contrast diagram 
on the (logNni, &)-plane in Figure |5| Note that, in this case, it is 
not the logarithmic, but the linear difference which is shown. The 
highest precision is of the order of 10"'* or even better. However, 
for the sake of simplicity, any value below 10"'' has been coded as 
zero. The largest discrepancy between both methods, i.e. the lower 
limit in the precision of our method if one takes Monaghan's as 
fiducial, amounts to 0.01, in agreement with the result shown in 
Figure El 

As can clearly be seen, the (lower limit in the) precision of 
our method depends both on b and Nhi. While the dependence on 
b extends to the whole range 10.0 - 100.0 km s"', the dependence 
on Nhi is limited to the range log Nhi = 16-20 dex. For a fixed b, 
there is a regime of values around 0.01 with a width of nearly 2 dex. 
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12 13 14 15 16 17 18 19 20 21 
log (Nhi) 



0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 

Figure 5. Worst Scenario: Lower limit to the precision of our method to 
synthesise Voigt profiles, given as the maximum value of the difference SV 
(cf. text for definition) between a Lya absorption profile computed accord- 
ing to equation Q using our approximation to the Voigt-Hjerting function 
(ea. l24l and l25l and Monaghan's algorithm, for the whole range of values 
for the parameters (log Nhi, fc) characteristic to intergalactic HI. The value 
con'esponding to each pair (log Nhi, is the maximum value of the quantity 
SV in the entire wavelength range A e [100, 2300] A. For clarity, any value 
below lO"'* has been coded as zero. 



: ^ 
: L 


r : 
J : 







1210 1212 1214 1216 1218 1220 



1210 1212 1214 1216 1218 1220 



r 



1200 1205 1210 1215 ^ 1220 1225 1230 

Wavelength [A] 

Figure 6. Examples of Worst Scenarios: Lyn absorption line profile synthe- 
sised using our method (upper panel in each row) and corresponding loga- 
rithmic difference with respect to our fiducial function (lower panel in each 
row) for the parameter pairs (logNnL^)) = (17.0, 10.0), (18.0,20.0), and 
(19.0, 70.0), from top to bottom, for which the maximum discrepancy with 
respect to the fiducial value of H in Figure|5]is largest. The maximum dif- 
ference of the order of 0.01 shown in Fi gurelslamounts to an extremely nar- 
row wavelength range between Gaussian core and Lorentzian wings. Away 
from these ranges, the accuracy improves significantly. For reference, and 
for comparison with Figurelsl a line corresponding to a constant logaiithmic 
dift'erence of -4 dex has been included in each of the lower panels. 

which gives as a result a narrow region of low accuracy all across 
the plane. Along this stripe, the difference reaches its highest value 
of the order of 0.01 or smaller for the combinations (low b, low 
Nhi) or (high b, high Nhi), in the b— and Nni-ranges stated above. 
Outside this stripe, the accuracy increases dramatically to values of 
the order of 10"'' or even better. 



The origin of any inaccuracy in our method is obviously the 
fact that terms of order n > 2 have been neglected in the series ^3), 
and furthermore, that the second term is this series has been ap- 
proximated as well. In particular, the origin of the 'low-accuracy' 
stripe on the (logNHi,i')-plane can qualitatively be understood in 
terms of the functional dependence of SV on b and Nhi- Consid- 
ering that both our method and Monaghan's take the zeroth order 
term exactly into account, it is legitime to state that the quantity 
A/f = H„^,y — Huon has a least a dependence of first order on a, i.e. 
AH(a, x) = a ■ h(a, x) where A is a function that may be of zeroth 
order in a. Hence, using equation Q and this last expression, it 
follows that 

AT(A) = C:~-h(a,x), with C: = —^ A^fiT,. 

As can be seen, t strongly depends on the ratio ^ . Therefore, an 
increase in b of one order of magnitude (from 10 to 100 kms"') 
is nearly compensated (in the sense that the value of r remains 
nearly constant) by an increase in Nhi of two orders of magni- 
tude, which accounts for the stripe of 2 dex in column density 
seen on the (log Nhi, fc)-plane. Why this happens precisely between 
log Nhi = 16-20 dex, as well as the shape of this stripe, is non- 
trivially related to the exact value of the constant C'., the fact that 
h may depend also on b through the damping parameter a, and the 
fact that the AV depends effectively not on At, but on e"'^'". 

It is worth mentioning that, for higher-order Lyman transi- 
tions, the precision of our method to compute Voigt profile has 
the same behaviour on the (logNHi,Z')-plane, and is the same as 
or even better than the precision of the Lyo- transition shown here. 
The reasons for this are, first of all, that the functional form of 
the absorption coefficient is obviously the same, irrespective of the 
transition. In addition, the lowest precision possible of 0.01 is the 
same for the whole range in a spanned by the Lyman transitions, 
according to Figure |4| Furthermore, higher transitions have lower 
damping parameters and our approximation is better the lower a, 
as already mentioned. Finally, since the constant C- is smaller the 
higher the order of the transition, the critical range of lowest pre- 
cision is shifted to higher column densities and higher Doppler pa- 
rameters. Since the ranges in Nhi and b are fixed for intergalactic 
HI, this has the net effect of increasing the high-precision region 
on the (log Nhi, b)-pla.ne-i.e. the region to the left of the stripe-for 
higher order transitions. In other words, the precision of our method 
improves from Lyo- to higher Lyman transitions. 

Even though a discrepancy of the order of 0.01 when using our 
method to model Voigt profiles may seem large, we want to empha- 
sise again that this is merely an lower limit for the precision in the 
entire wavelength range A e [100,2300] A, in the case of Lyo-. 
It turns out that the range in wavelength for which the accuracy is 
lowest is negligible for practical purposes. In order to show this, we 
first choose three points on the (log Nhi, ii)-plane along the stripe of 
lowest precision, i.e. for which the maximum difference is largest. 
Using these parameters, we synthesise Lycf absorption profiles us- 
ing our method and Monaghan's, and compute again for each of 
these "worst scenarios" the quantity SV for the whole wavelength 
range. The result is shown in Figure|6| Each row corresponds, from 
top to bottom, to the parameter pairs (logNHi,^) = (17.0, 10.0), 
(18.0, 20.0), and (19.0, 70.0) chosen along the low-precision stripe. 
The upper panel of each row shows the hya absorption profile syn- 
thesised using our method, and the corresponding lower panel gives 
the logarithmic difference log,Q SV between our method and Mon- 
aghan's as a function of wavelength. As can be seen, the small- 
est discrepancies are given at the line cores, as expected, since in 
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this regime the zeroth order term dominates and both our approach 
and Monaghan's exactly take this term into account. The largest 
discrepancies, of the order of 0.01, are present in an extremely 
narrow range of ^ 0.06 A for (logNHi,/?) = (17.0, 10.0), of 
AA « 0.17 A for (logNm,/?) = (18.0,20.0), and of A^ « 0.53 A 
for (logNHi,^) = (17.0,20.0). This discrepancies are found at the 
boundaries between Gaussian core and Lorentzian wings, due to 
the fact that our method neglects terms of order n > 2, which dom- 
inate the behaviour of H in that regime. Note, however, that the 
difference rapidly drops with increasing distance (in A) from the 
line center to values of the order of e.g. 10"'" at a distance A/1 « 5 
A for the first two rows, and AA x 15 A for last row. Hence, the 
effective accuracy of our method is far better than 0.01 in the wave- 
length range shown here. For the same reason mentioned in the last 
paragraph, the precision for higher order Lyman transitions is of 
the same order as or even better than for the Lyo- transition shown 
here. 



6 APPLICATION 

As a further test of the quality of our approximation to model Voigt 
profiles, and to show its accuracy in a less academic situation as 
in the last section, we consider fitting a synthetic spectrum to a 
real quasar absorption spectrum with a population of intergalactic 
HI absorbers spanning a representative range in h and Nhi along a 
random line-of-sight. For this purpose we use the observed spec- 
trum of the quasistellar source QSO J2233-606, a relatively bright 
(B = 17.5) quasar at an intermediate redshift Ze,,, = 2.238. 

The spectrum of the source QSO J2233-606, centered at the 
HDF-S, was obtained during the Commissioning of the UVES in- 
strument at the VLT Kueyen Telescope and reduced at the Space 
Telescope European Coordinating Facility. The high-resolution 
spectroscopy (R « 45000) was carried out with the VLT UV- 
Visual Echelle Spectrograph (UVES). The data were reduced in the 
ECHELLEAJVES context available in MIDAS. The final combined 
spectrum has constant pixel size of 0.05 A and covers the wave- 
length range 3050-10000 A. The S/N ratio of the final spectrum is 
about 50 per resolution element at 4000 A, 90 at 5000 A, 80 at 6000 
A, 40 at 8000 A. The data used here are publicly available and were 
retrieved from www.stecf.org/hstprogrammes/J22/J22.html, in its 
version of November 23, 2005. 

With help of th e MIDAS package FITLYMAN 
jpontana & BallesteT 1995), which performs line fitting through 
X'^ minimization of Voigt profiles, Cristiani & D'O dorico ( 2000) 
determined the redshifts, column densities, and Doppler widths of 
the identified absorption features imprinted in the spectra of QSO 
J2233-606. In this way, they found that the line of sight to QSO 
J2233-606 intersects a total of 270 Lya Forest clouds, and identi- 
fied other 24 absorption systems containing metal lines. The Lyo- 
absorbers span a range in column density of 10'-^" - 10" '" cm"^, 
and a range in Doppler parameters of 1.0 - lll.Okms"'. The 
distribution of the parameter pairs (Nhi,^) for these systems for 
be [10, 100] km s"' is shown in Figure^ 

Using this list of hya line parameters (Zats, b, Nhi), we gener- 
ate a synthetic spectrum of the QSO J2233-606 in the wavelength 
range 3340 - 3530 A with a higher resolution than that of the ob- 
served spectrum of AA = 0.01 A, using eqs. Q, <23> . <24> . and 
<25> . We synthesise a second spectrum using Monaghan's algo- 
rithm, and compute the logarithmic difference 6V between these 
two synthetic spectra in the same fashion as in the previous sec- 
tion. In this way, we test again our method against the highest- 
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Figure 7. Distribution in column density and Doppter parameter of the 
absorbing systems along the line-of-sight towards the source QSO J2233- 
606. This particular LOS contains a total of 264 absorbing systems with 
b £ [10, 100] kms"' in the range 3050 - 10000 A (open circles), 69 from 
which are in the wavelength range 3340 - 3530 A (filled circles). The data 
was taken from Cristiani & D'Odorico ( 2000). For reference and compari- 
son, the contrast diagram shown in Figure|5]has also been included. 

precision method available, for a typical range in column densities 
and Doppler parameters present in a QSO spectrum. The result is 
shown in Figure[8| The upper panels of each row show the observed 
spectrum and the spectrum synthesised using our method, whereas 
the lower panels show the logarithmic difference between both syn- 
thetic spectra. We choose to cut off the logarithmic difference at 
10", since differences smaller than these are out of the range of 
the highest available numerical precision. It can be seen again, as 
in Figure |6| that the smallest discrepancies are given at the line 
cores, and the largest, of the order of 10"'', are given at the wings 
(cf. discussion of Figure |S| Section l53l . As can be seen from the 
column density and Doppler parameter distribution in FigureQ the 
largest discrepancies in this wavelength range are consistent with 
the maximum absolute differences shown in Figure |5| Note that 
in the spectral regions where no apparent absorption features are 
found, the logarithmic difference does not fall to -oo, as one would 
naively expect. These features are present in Figure|6|as well. The 
reason for this 'valley-shaped' features is that, even though having 
small values away from the line center, the function H does not fall 
to zero, and thus in these regions the wings of two or more lines 
overlap. Strictly speaking, in these regions there is always some 
absorption left, i.e. t < 1, and different methods to compute H will 
account for this effect differently. Again, since our method neglects 
terms of order n > 2, the absorption in these regimes differs from 
its fiducial value. The lack of this terms in our approximation to H 
is also pointed out pictorially by the local maxima symmetrically 
placed around the deeps corresponding to the logarithmic differ- 
ence at the line cores. 



7 SUMMARY 

The absorption lines imprinted in the spectra of background sources 
yield a wealth of information about the physical and chemical prop- 
erties of the intervening absorbing material, as is the case of inter- 
vening neutral hydrogen (HI) systems embedded in the intergalac- 
tic medium (IGM) and associated with galaxies and larger struc- 
tures. In order to extract the desired information from these absorp- 
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Figure 8. Observed (dashed line) and synthetic (solid line) spectrum of the quasar HDF-S QSO J2233-606 (upper panel in each row) in the wavelength 
range [3340, 3530] A ( cf. text for reference). The synt hetic spectrum was generated using our approximation to the Voigt-Hjerting function and the list of line 
parameters obtained bv lCristiani & D'Odoricol <200Cl) . The lower panels show the logarithmic difference between our synthetic spectrum and one generated 
using the same list of line parameters and Monaghan's algorithm to calculate Voigt profiles. For reference, we include in the difference panels an horizontal 
line conesponding to a logarithmic difference of -4 dex. 



tion lines, their profiles have to be modeled in a proper way. In the 
case of absorption features found on QSO spectra, absorption line 
profiles are best modeled by Voigt profiles, which are mathemati- 
cally given in terms of the Voigt-Hjerting function. 

In this work, we presented a simple analytic approximation to 
the Voigt-Hjerting function with which Voigt profiles can be mod- 
eled for an arbitrary range in wavelength (or frequency), column 
densities up to 10"^ cm"^, and for damping parameters satisfying 
a < lO""*. Starting with an exact expression for the Voigt-Hjerting 
function in terms of Harris' expansion that is valid for a < 1, we 
showed that the zeroth order term of this series, a Gaussian func- 
tion, is suitable for modeling absorption line profiles emerging in 
a medium with low column density Nhi g 10"cm"^. However, for 



higher column densities, terms of higher order have to be taken 
into account. A key point leading to this conclusion is the fact that 
it is not the damping parameter alone, but rather the factor a ■ Nhi 
that determines to which extent terms of order higher than zeroth 
in Harris' expansion may or may not be neglected. We showed that 
the departure of the actual Voigt-Hjerting function from the first 
two terms in Harris' expansion is of the order of lO"' or less for 
an arbitrary wavelength range and a < 10"''. Hence, we concluded 
that with an approximation to first order in a to the Voigt-Hjerting 
function Voigt profiles can be modeled with moderate to high ac- 
curacy. 

On this basis, we obtained a simple analytic expression for the 
Voigt-Hjerting function and consequently for the absorption coeffi- 
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cient of intergalactic HI, in terms of an approximation for tlie sec- 
ond term Hi of Harris' expansion. The main advantages of tlie ana- 
lytic expression we presented here are, first, that it is valid for an ar- 
bitrary wavelength range, in the sense that is has no singularities. In 
addition, it is simple and flexible to implement in a variety of pro- 
gramming languages to numerically compute Voigt profiles with 
moderate speed and moderate to high accuracy. As a matter of fact, 
our method to compute the Voigt-Hjerting function is faster with re- 
spect to other known standard methods, for instance, Humlfcek's or 
Monaghan's algorithm. Furthermore, our approximation reaches an 
accuracy of 10"'' or better in a wide wavelength range, and of the 
order of than 10"^ only a negligible wavelength interval, for val- 
ues of a and Nhi characteristic to intergalactic HI absorbers. Our 
method thus offers a great compromise between speed, accuracy, 
and flexibility in its implementation. 

Even though we did not extend our discussion in this work 
to other transitions typically present in quasar absorption spectra 
and associated to HI absorbers, such as metal lines, our method to 
synthesise Voigt profiles can certainly be applied to most of these 
elements as well, since their column densities are obviously the 
same, and their ranges in a strongly overlap with the range of a for 
intergalactic HI, for which our approximation to the Voigt-Hjerting 
function is valid. As a matter of fact, our approximation is valid 
to model absorption Voigt profiles found in any type of spectrum 
(stellar, solar, etc.), which arise in a medium whose damping pa- 
rameter and column density satisfy a < lO""* and Nhi < 10"^ cm"^, 
respectively. 
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